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1 Introduction 

This work has been supported in part by the Netherlands Organization for Scientific Research (APERTIF 
project funded by NWOGroot), the Swedish Agency for Innovation Systems VINNOVA and Chahiiers Uni- 
versity of Technology (VINNMER - Marie Curie Actions international qualification fellowship). Some results 
of this work have been recently published in [3 [2] . 

This report is a shortened version of the original ASTRON report (RP324, 2011) which is publically 
availible upon request. The purpose of this report is to document the noise performance of the Apertif 
array receiver which has been used to measure its performance as an aperture array and to characterize 
the recently developed noise measurement facility THACO. The receiver system is the second prototype of 
APERTIF (DIGESTIF2) [3J[1] that includes the array antenna of 144 dual-polarized TSA elements, 144 Low 
Noise Amplifiers (LNAs) (Tmin =35-40K) and the data recording/storing facilities of the initial test station 
that allow for off-line digital beamforming. 

The primary goal of this study is to compare the measured receiver noise temperatures with the simulated 
values for several practical beamformers, and to predict the associated receiver noise coupling contribution, 
antenna thermal noise and ground noise pick-up (due to the back radiation). The measurements were 
performed over a wide frequency band and scan range. In the course of the measurements, the 4, 16, 25 
and 49 active antenna elements were used in beamforming, while the remaining elements were connected to 
the LNAs with the outputs terminated in matched loads. The experimental results were obtained when the 
receiver with 4- and 16-element analog beamformers was placed inside THACO (a metal shielding cabin which 
was designed to mitigate the noise contributions due to the ground and trees) and in the open environment 
at Westerbork near the Westerbork Synthesis Radio Telescope (WSRT), using digital beamforming. The 
details of these measurements can be found in ^ . 

Since, the digital data recording-storing facilities are available for the DIGESTIF array receiver, one 
can also evaluate the system noise temperature for beamforming scenarios that are based on standard signal 
processing algorithms [S]-[B] and the noise correlation cocfBcicnts between the array channels. In this report, 
we therefore show the range of the realized system noise temperatures when an embedded element only and all 
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array elements are used in beamforming. These results demonstrate the pros and cons of the considered low- 
gain antenna and high gain digital array receivers for the purpose of the noise temperature characterization 
in an open environment and inside a shielding cabin such as THACO. 



2 The model of the array antenna 

The array consists of 2 x 72 aluminium Tapered Slot Antenna (TSA) elements with a pitch of 11 cm (0.52A 
at 1420 MHz) on a rectangular 8x9 grid. Each TSA is fed by a wideband microstrip feed which has been 
integrated with an LNA on a printed circuit board. This design features a very short transmission line 
between the antenna and the LNA, since the circular slotline cavity has been moved sideways. More details 
on the array design and the numerical approach used for the EM-analysis of the antenna can be found in [7] 
and [6l [8] , respectively. The simulations of the antenna have been carried out using CAESAR software that 
is an array system simulator, developed at ASTRON [^. Note that these simulations did not take into 
account the effect of the contributions from obstacles (trees and telescopes/buildings) near the horizon. 

3 The noise model of the receiver system 

The noise model used in this study is based on the equivalent system representation as described in [101 . 
According to this representation, the sensitivity of the array receiver can be computed as follows: 



where and Tgys are the physical antenna area and system noise temperature, respectively. The latter 
consists of three main contributions: (i) The external noise contribution - the ground noise picked up due to 
antenna back radiation, which was computed from the simulated illumination pattern of the antenna array 
for the specified beam former weights, (ii) the thermal antenna noise due to the losses in the conductor 
and dielectric materials of TSAs and microstrip feeds. The conductor losses are computed through the 
evaluation of the antenna radiation efficiency using the methodology detailed in [11] and the dielectric losses 
are computed based on the experimental evaluation of the feed loss [H] and (iii) the noise due to LNAs 
which is dependent upon the noise properties of LNAs and active reflection coeflticients seen at the ports 
of the antenna array. Note, that by using this definition of Tsys, all noise temperature contributions are 
referenced to the sky (in front of the antenna aperture). 

It is important to note that the system noise temperature which is calculated from the measured Y- factor 
is also referenced to the sky. Therefore, to distinguish between the external noise, antenna thermal noise and 
receiver contributions, one needs to know the beam shape (to compute Text) and radiation efficiency jyrad of 
the antenna, which are in general dependent on the beamformer weights. 

4 The model of the beamformer used to compute the optimal 
weights 

The beamformer model used in this study is based on the mathematical framework which has been described 
in [6]. Within this framework, the array receiver system is subdivided into two blocks: (i) the front- 
end including the array antenna. Low Noise Amplifiers (LNAs); and (ii) the beamformer with complex 
conjugated weights {w*}^^]^ and an ideal (noiseless/reflectionless) power combiner realized in software. 
Here, = [wl, . . . , w^] is the beamformer weight vector, H is the Hermitian transpose, and the asterisk 
denotes the complex conjugate. Furthermore, a — [ai, . . . , a^]'^ is the vector holding the transmission-line 
voltage- wave amplitudes at the beamformer input (the N LNAs outputs). Hence, the fictitious beamformer 
output voltage v (across Zq) can be written as w = w^a, and the receiver output power as |wp ~ vv* = 




(1) 
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(w^a)(w'^a)* ~ (w^a)(a"^w*)* — w^aa^w, where the proportionahty constant has been dropped as this 
is customary in array signal processing and because we wiU consider only ratios of powers. 

Although each subsystem can be rather complex and contains multiple internal signal/noise sources, it 
is characterized externally (at its accessible ports) by a scattering matrix in conjunction with a noise- and 
signal- wave correlation matrix. In this manner, the system analysis and weight optimization becomes a purely 
linear microwave circuit problem. The sensitivity metric ^cff/7sys, which is the effective area of the antenna 
system divided by the system equivalent noise temperature, can be expressed in terms of the Signal-to-Noise 
Ratio (SNR) and the normalized flux density S'somce of the source (in Jansky, 1 [Jy] = IQ-^e [Wm-^Hz"^]) 
as 

^'^^ ''^-SNR, where SNR = 1?^, (2) 



>- sys ^source 



and where fc^ is Boltzmann's constant. The SNR function is defined as a ratio of quadratic forms where C 
is a Hermitian spectral noise-wave correlation matrix holding the correlation coefficients between the array 
receiver channels, i.e., Cmk = E{cmC^.} = Cmc\. (for A;,m = 1 . . .N). Here, Cm is the complex-valued voltage 
amplitude of the noise wave emanating from channel m (see [13] and references therein) , which includes the 
external and internal noise contributions inside the frontend block. We consider only a narrow frequency 
band, and assume that the statistical noise sources are (wide-sense) stationary random processes which 
exhibit ergodicity, so that the statistical expectation can be replaced by a time average (as also exploited 
in hardware correlators). C is nonzero if noise sources are present in the external environment and inside 
the system, due to e.g. the ground, LNAs, and sky. For a single point source on the sky, the signal-wave 
correlation matrix P = ee^ is a one-rank positive semidefinite matrix. The vector e = [ei, 62, .., cat]'^ 
holds the signal-wave amplitudes at the receiver outputs and arises due to an externally applied plane 
electromagnetic wave E\. 



4.1 Maximum Sensitivity 

Maximizing Eq. ([2]) amounts to solving the largest root of the determinantal equation [5]: det (P — SNRC) = 
(c/. |14j). Next, the optimum beamfomer weight vector WMaxSNR is found through solving the correspond- 
ing generalized eigenvalue equation P WMaxSNR = SNR C WMaxSNR for the largest eigenvalue (SNR) as de- 
termined in the previous step. The well-known closed-form solution for the point source case, where P is of 
rank 1, is given by |14) 

WMaxSNR = C~^e, with SNR = e-'^WMaxSNR (3) 

where the eigenvector e corresponds to the largest eigenvalue of P. 



4.2 Maximum output power or the Conjugate Field Matching (CFM) condition 

When C equals the identity matrix I (thus equal and uncorrelated output noise powers), the receiver output 
noise power w^Cw = w^w becomes independent of w in case its 2-norm w^w is a constant value, typically 
chosen to be unity. With reference to ([3]), the weight vector that maximizes the received power, and thus 
realizes a maximum directive gain (and effective area) in the direction of observation, is therefore 

wcFM = e. (4) 

These weights optimally satisfy the Conjugate Field Matching (CFM) condition [HI [HI [17] . 



4.3 Minimum system noise temperature 

Similarly, one can develop an expression for computing the beamformer weights for the minimum Tsys , that 
is the case when the source of interest has no contribution and thus independent on the weights. For this 
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case, the optimal beamformer is described as: 



WMinTsys — (5) 

where Cn = 1. 



5 Numerical results for the 144-channel (full-polarization) beam- 
former 

5.1 Simulation details 

The simulations were performed with the newly developed numerical tool box for the CAESAR software 
[5]. This toolbox was initially aimed at the analysis and optimization of the PAF systems, and has been 
interfaced with GRASP to compute the overall noise wave scattering matrix due to external and internal 
noise sources as well as the secondary array patterns after the scattering from the dish. For the present 
study, we have used the pre-processor of this software to determine the optimal beamformer weights (so 
as to account for the non-uniform noise distribution of the environment) and to evaluate the receiver noise 
contributions due to internal noise sources according to the model presented in [TD] and 
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5.2 The array beam noise temperature and its contributions 
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Figure 1: (a) The simulated array beam noise temperature versus frequency and its contributions due to (b) 
the noise couphng effects in the receiver, (c) the ohmic losses of the antenna and (d) the external noise due 
to the ground (due to the back radiation of the array, when Q > 90°). 



6 Results for experimental (bi-scalar) beamformers with 4, 16, 25 
and 49 channels 

This chapter shows the simulated array antenna properties and compares them with measured receiver 
noise temperatures versus frequency and scan angle for various beamformer configurations. It starts with 
a description of the measurement setups and method. It then presents the simulation results for the array 
antenna patterns and the array beam noise temperatures and its contributions. Finally the results of cross- 
comparison of the measured Tree at Westerbork and inside THACO are presented. 
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6.1 Measurement method and setups at Westerbork and inside THACO 

For the noise measurements of the APERTIF tile as an aperture array the Y- factor hot / cold method has been 
used |18]-|19|. During the measurement the tile is placed horizontally on the ground for the measurement in 
Westerbork (in the open area) or inside the big shielding facility for the Y- factor hot/cold measurements that 
is called THACO. Two methods were followed, one using analog beam forming with 2x2 and 4x4 elements for 
the measurements inside THACO, the other using digital beam forming for various beam configurations at 
the WSRT location. The output signals from the LNAs of a 2x2 and 4x4 array in the centre of the tile inside 
THACO were added with in-phase combiners, forming broadside beams. The analog output signals of the 
beam formers were fed to the input of an Agilent Noise Figure Meter 8970B and the noise temperature was 
determined with the Y-factor method, using the cold sky as a 'cold' load and the roof of THACO, covered 
with absorbing material, as the 'hot' load. The digital beam forming and processing of the APERTIF 
prototype system at the WSRT provided a much more flexible system, with which beams could be formed 
with a larger number of elements, pointing in any desired direction. For the digital processing method, a 
total of 49 individual antenna elements and LNAs (limited by the number of available receivers at the time 
of measurements) are connected via 25 m long coaxial cables to the back-end. The back-end electronics 
is located in a shielded cabin, together with the down converter modules and digital processing hardware 
[20] , Data are taken with the array facing the (cold) sky as 'cold' load, after which a room temperature 
absorber is placed over the array for the measurement with the 'hot' load. The data processing takes into 
account correlations between data from individual elements. The results are stored in a covariance matrix as 
a function of frequency. Using off-line digital processing, beams with a combination of any of the 49 active 
elements can be formed and beams may be scanned in any direction by applying weights to the elements of 
the covariance matrix. In this way the equivalent beam noise temperature as a function of frequency from 
1.0 to 1.8 GHz has been determined for 2x2, 4x4, 5x5 and 7x7 element arrays, looking at broadside. Also 
the equivalent beam noise temperatures as a function of scan angle for the 4x4 and 7x7 element arrays have 
been determined. 

6.2 Simulated antenna beam directivity and noise temperatures 

In this section, we present the numerical results for practical beamformers for the frequencies ranging from 
1 GHz to 1.6 GHz and over the scan range within which 9 changes from 0° to 85° and = — 360°. 
The practical beamformers combine the signals received by 4, 16, 25 and 49 elements. The numerical 
results include the directivity of the array antenna (see Fig. [2] and the receiver noise temperature and its 
contributions due to noise coupling effect, antenna ohmic loss, and ground noise due to back radiation (see 
Fig.!-??). 

Figure Ela) shows the receiver noise temperatures of the DIGESTIF tile that were computed for the 
boresight direction of observation at 8 frequency points within the bandwidth of 1-1.6 GHz. These results 
clearly demonstrate that for all practical beamformers the on-axis beam noise temperature is weakly depen- 
dent on frequency and takes values between 42 and 61 K that are 5-30% higher than the minimum noise 
temperature of LNAs (Tmin = 35 — 40)K. The noise contributions due to the receiver noise coupling effects 
Tcoup, antenna ohmic losses Tiad and external (ground) noise pick-up Text, as shown on fig. Elb), ^c) and 
ISjd) do not exceed 13, 3.5 and 6.5 K, respectively. At 1.42 GHz - the frequency at which the array design was 
optimized - the temperatures Tcoup and Tf^xt and the total receiver temperature take the minimum values 
within the operational bandwidth as the result of the minimized impedance mismatch loss and relatively 
low side and back radiation levels. Figure |4lja)-(d) shows how the receiver noise and its weight-dependent 
noise components vary with scan angle for beamformers with 16 and 49 elements at 1 GHz and 1.42 GHz. 
For these beamformers, respectively, the increase of the noise temperatures is less than 20% when the scan 
angle is smaller than ~ 30° and ^ 40° off boresight direction. For larger scan angles, however. Tree rapidly 
increases and becomes as high as 80-160 K depending on the number of active antenna elements, scan plane 
and frequency. Such high values are mainly due to the strong mutual coupling between antenna elements 
at low frequencies (causing the rise of the receiver noise coupling contribution as observed on fig. [SJc)) and 
high side-lobe levels at high frequencies for scanned beams. 
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Figure 2: The array directivity versus frequency. 



6.3 Comparison with measurements in the open environment at Westerbork 

This subsection compares the predicted and measured receiver noise temperatures versus frequency. The 
results of cross-comparison of the measured Tree at Westerbork (in the open environment) and inside THACO 
(which is expected to shield the receiver from the ground noise) are presented. Also, the measured and 
modeled noise temperatures of a single Vivaldi element receiver inside THACO are shown. 

Upon collecting the simulation and measurement results, we can compute the relative difference between 
the modeled and measured receiver noise temperatures. Figures [6][7] show both the simulated and measured 
Tree as well as their relative difference as a function of frequency and scan angle. As one can see on figlH^c), 
the relative difference between simulations and measurements is smaller than 20% over the entire frequency 
range for all practical beamformers, except for the beamformer with 4 active channels in the region of 1-1.2 
GHz. At these frequencies the 4-element subarray has a rather low gain (<10 dB) and, thus the experimental 
receiver picks up the noise due to the buildings and trees that are present in the actual environment, but 
were not accounted for in the model. This reasoning is supported by the measurements which were done 
inside the shielding cabin (THACO) (see figlUJc)). For the latter tests, the measured noise temperatures 
for the 4-channel beamformer lie within the region of the simulated values with 15-20% difference at most 
frequencies. 

The agreement between the simulated and measured noise temperatures over the scan range is also good, 
and for the beamformers with 16 and 49 channels was found to be at the level of <25% within the scan range 
of ±(40 — 45°). For larger off-axis angles, the measured temperatures are much higher with the difference 
up to 50% relatively to the corresponding simulated values. 
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Figure 3: (a) The simulated receiver noise temperature versus frequency and its contributions due to (b) 
the noise couphng effects in the receiver, (c) the ohmic losses of the antenna and (d) the external noise due 
to the ground (due to the back radiation of the array, when 6 > 90°). All temperatures are for broad side 
beams. 



6.4 Comparison with measurements inside THACO 
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Figure 4: The simulated receiver noise temperatures versus scan angle in three scan planes for beamformers 
with (a,b) 16 and (c,d) 49 active antenna elements at 1.0 GHz and 1.42 GHz. 



7 Conclusions 

The noise performance of the aperture array tile receiver has been simulated and compared to the experi- 
mental results as obtained though the 'hot-cold' measurement procedure inside a shielding cabin (THACO) 
and in the open environment near WSRT. The measurements have been carried out for several practical 
beamformers (with 4, 16, 25 and 49 active channels) over the frequency band from 1 to 1.6 GHz and a 
wide 3D beam scan range. The presented numerical results include the antenna patterns and system noise 
temperatures Tsys for all considered practical situations, as well as separate noise contributions due to the 
receiver noise coupling effects Tcoup, antenna ohmic loss Tiad and external (sky and ground) noise Text- 

The numerical results demonstrate that the on-axis beam noise temperatures take values ranging be- 
tween 42 and 61 K within the frequency bandwidth and are maximum 30% higher than the minimum noise 
temperature of LNAs {T-adn = 35 — 40)K. The noise contributions Tcoup, ?rad and Text do not exceed 13, 
3.5 and 6.5 K, respectively for broad side beams. At 1.42 GHz - the frequency at which the antenna was 
optimized - Tsys is lowest, as the result of the minimized impedance mismatch loss and relatively low side 
and back radiation levels. 
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(b) (c) 

Figure 5: (a) The simulated receiver noise temperature ([K]) over scan angle at 1.0 GHz and its contributions 
due to (b) the noise couphng effects in the receiver, (c) the ohmic losses of the antenna and (d) the external 
noise due to the ground (due to the back radiation of the array, when 9 > 90°). 



The receiver noise temperature exhibits a strong dependence on the beamformer weights and degrades 
when scanning far off boresight direction. For beamformers with 16 and 49 channels respectively, the relative 
increase of Tree was found to be less than 20% for the scan angles smaller than ^ 30° and ~ 40°, and a 
factor 2-4 for larger angles, depending on frequency and beamformer. Such high values of Tree are mainly 
due to the strong mutual coupling between antenna elements at low frequencies (causing the rise of Tcoup 
and high side-lobe levels at high frequencies for scanned beams. 

There is a good agreement between simulations and measurements that were performed in the open 
environment at Westerbork. The relative difference between the modeled and measured Tsys is smaller 
than 20-25% over the entire frequency band and within the scan range of ±(40 — 45°), except for the 4- 
channel beamformer for which this difference can be twice as large. For the latter beamformer case, the 
antenna pattern is rather broad (the gain is <10 dB) and, thus the experimental receiver can pick up an 
additional noise component due to the buildings and trees that are present in the actual environment, but 
were not accounted for in the model. This reasoning is supported by the measurements inside the shielding 
cabin (THACO) for which the agreement with simulations significantly improves. Furthermore, for large 
off boresight scan angles, the temperatures as measured in the open environment are much higher than the 
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Figure 6: (a) The simulated and (b) measured system noise temperatures (Tsys = T^-cc + ^Isky) versus 
frequency and (c) the relative difference between them. 



simulated ones, most likely due to the above mentioned effects of the noisy environment of which the exact 
temperature distribution is not well known. 
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